Ultrafast (but Many-Body) Relaxation in a Low-Density Electron Glass 
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We present a study of the relaxation dynamics of the photoexcited conductivity of the impurity 
states in the low-density electronic glass, phosphorous-doped silicon Si;P. Using subband gap opti- 
cal pump-terahertz probe spectroscopy we find strongly temperature and fluence dependent glassy 
power-law relaxation occurring over sub-nanosecond time scales. Such behavior is in contrast to the 
much longer time scales found in higher electron density glassy systems. We also find evidence for 
both multi-particle relaxation mechanisms and/ or coupling to electronic collective modes and a low 
temperature quantum relaxational regime. 
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Glasses are ubiquitous states of matter with position- 
ally or rotationally randomized frozen degrees of free- 
dom. Their inherently frustrated interactions give rise to 
the existence of many energetically lovi^-lying metastable 
states, which have a wide distribution of potential barri- 
ers separating them in configuration space. At low tem- 
peratures, such systems are kinematically constrained 
from reaching their true ground-state on typical experi- 
mental time scales. They are characterized by long re- 
laxation times, memory, and aging effect^^"'^. 

In analogy with structural glasses, electronic glasses 
can be defined as systems with a random distribution 
of localized charges^^-. Here long relaxation times and 
glassy phenomena derive from a combination of dis- 
order and long-range unscreened Coulomb interaction. 
The motion of any one charge manifestly necessitates a 
many-particle rearrangement of the other electron occu- 
pations to reach lower energy. Such many-particle pro- 
cesses are inherently slow and inefficient at finding the 
true ground-state configuration. Electronic glasses or 
'Coulomb glasses' may be realized in granular metals and 
amorphous and lightly doped semiconductors, which are 
all expected to exhibit certain similar qualitative behav- 
ior. In addition to explicit glassy effects, there are also 
predictions for their equilibrium properties like the DC 
hopping conductivit ySI and the power-law dependencies 
of the AC response2ni21. 

Recently there has been great progress in the under- 
standing of how glassy non-equlibirum effects manifest 
themselves in these systems. In amorphous compounds 
like In2 03_2: the natural history- free relaxation law of 
the DC conductivity after excitation is logarithmic over 
more than five decades in time^, a behavior that may 
be rooted in a broad distribution of relaxation times. In- 
teresting temperature dependencies have been found at 
low temperature, which have been ascribed to qua ntum 
tunneling instead of thermally activated hoppingl^. 

The majority of such experiments have been on 
syste ms such as microcrystalline In2 03_2; or granular 
metala^E^E^in which the charge density is high. It would 
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be interesting to look for such phenomena in systems 
with much lower densities like doped semiconductors. Al- 
though the equilibrium transport properties expected for 
electro nic glasses have been reported in systems like Si:B 
Qj. g^.p| io | ii | i6 | i7 |^ g^g fg^j. g^g .^g ]jjjQ-^y^ evidence for true 

glassy behavior in the conductivity, i.e. long-time non- 
exponential relaxation indicative of metastability and 
frustration has not been observed. This is likely due 
to their much faster time scales owing to the much 
smaller charge densities, which is crucial for their quan- 
tum dynamics^^. So while in principle glassy relaxation 
should exist in such systems, it is expected to be too fast 
to be resolved with usual techniques. 

In this Letter we report the use of ultrafast optical 
pump-terahertz probe (OPTP) spectroscopy to resolve 
glassy relaxation in the doped semiconductor Si:P. Ter- 
ahertz time-domain spectroscopy is an ultrafast optical 
technique in which electric field transients are used to 
measure the conductivity of a material. Here a sample 
is optically excited with a laser pulse and then probed 
at a later time with a terahertz (THz) pulse to mea- 
sure the induced conductivity changes with picosecond 
resolution. While OPTP has been successfully applied 
to correlated electron syst ems such as high-Tc's, nian- 
ganites and semiconductoralSlini jn unraveling the various 
dynamics by their different relaxation time scales, this is 
the first such study on glasses. In this work, we find 
strongly temperature and fluence dependent glassy non- 
exponential relaxation occurring on a sub-nanosecond 
time scale, which is in contrast to the much longer time 
scales observed in high electron density glassy systems. 

Experiments were performed on nominally uncompen- 
sated phosphorous-doped silicon samples, which were cut 
from a Czochralski grown boule grown by Recticon En- 
terprises Inc. to a specification of 5 cm in diameter with a 
P-dopant gradient along the axis. This boule was subse- 
quently sliced and then polished down to 100 /im. Sam- 
ples from this boule have previously been used for an 
extensive study of the THz-range conductivity in the 
phononless regim^i^^i^l fn this study, we concentrated 



on a particular Si:P sample with a 300 K resistivity of 
0.0202 fi-cni. According to the Thurber scald^this cor- 
responds to an estimated phosphorus concentration n of 
1.37 X 10^^ dopant atoms/cm"^, which puts it at 39% of 
the way to the MIT. This sample's^ localization length 
(« 13 nm) is slightly larger than the interdopant spac- 
ing of w 9 nm and its static dielectric constant (w 14) 
is enhanced over that of undoped silicon (11.7) demon- 
strating that collective quantum effects may play a role. 
Although nominally uncompensated, it is also generally 
believed that this doping level is high enough that the 
material self-compensates i.e. the long-range Coulomb 
interaction ionizes a substantial fraction of P sites. 

The experiments utilized a regeneratively amplified 
Ti:Al2 03 laser system operating at 1 kHz, producing 
nominally 1.0 mJ, 150 fs pulses at 1.5 eV. The Si:P sam- 
ple was excited at 1500 nm (0.82 eV) via an optical para- 
metric amplifier with pump fluences ranging from 4.2 to 
424.4 ^J/cm^ producing photoexcited carrier densities^^ 
ripe ~ 3.2x10^^ — 3.2x10^^ cm^-^. The experiments were 
performed in transmission (T) with the Si:P sample in- 
side an optical He cryostat capable of reaching 4 K. The 
THz 'probe' pulses were generated and detected using 
electro-optic techniques. A schematic of the THz setup 
is given in Ref . 22 , Fig. [pi shows the transmitted electric 
field of the THz pulse before E{t) and after AE{t) opti- 
cal excitation. The decrease in the transmitted electric 
field is associated with the change in conductivity (i.e.. 
Act oc —AT/T). From an analysis of the skin-depth, the 
sample is presumed to be homogeneously excited at low 
temperature. In equilibrium, the THz range conductivity 
in these materials derives from transitions between local- 
ized levels and is believed to be proportional to a high 
power of the localization length ^ as cri ex iVg^^ where 
A^o is the non-interacting density of states^^. 

The excitation process is shown in Fig. [TJd. Electrons 
are excited from the impurity band to a relatively large 
energy in the conduction band (0.8 eV). However, they 
quickly decay back into impurity states d through succes- 
sive many-body hops find lower energy orbitals, slowly re- 
ducing the transmission. In order to probe only impurity 
band states it is essential that the pump laser's energy is 
below the band gap threshold (1.5 eV). Note that the dy- 
namics of conduction band-valence band recombination 
is completely different as shown in Fig. 111;. Recombina- 
tion of 1.5 eV electron- hole pairs produced by above-gap 
excitation occur on time scales of millisecond or longer, 
while 0.8 eV excitations decay on sub-nanosecond time 
scales. While recombination of conduction band - valence 
band pairs from a 1.5 eV excitation pulse is slow it can 
be described by a simple exponential, and is not glassy. 

In Fig. [2] we show the peak normalized photoinduced 
AT/T dynamics at various npe's at T = 5 K with 0.8 
eV excitation. All the data present a fast ^1 picosecond 
rise, a plateau, and then a much slower decay. As we 
will discuss below, the data have a number of unusual 
features consistent with glassy dynamics. At high npe 
the decay is power-law with an exponent of ^^-2 until 
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FIG. 1: (color) (a) Electric field E{t) of transmitted THz 
pulse, and induced change in electric field —AE{t) at T = 5 K. 
(b) Schematic energy diagram of Si:P with impurity states 
close to the conduction band. The upward arrow indicates 
photoexcitation at 0.82 eV. The downward arrows indicate 
the subsequent relaxation processes, (c) Comparison of the 
normalized differential transmission of the THz probe pulse 
peak AT/T versus delay with a pump excitation energy above 
(1.5 eV) and below (0.82 eV) the band gap at T = 5 K for a 
photoexcited carrier density of 3 x lO^'^ cm~^ in both cases. 
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FIG. 2: (color) Time-resolved AT/T photoinduced transmis- 
sion dynamics of Si:P sample at different npe's at T = 5 K 
and an excitation energy of 0.82 eV. The data are normalized 
to maximum signal. Fits are shown as black dashed fines. 
The inset shows the photoexcited charge density npe and flu- 
ence dependence of the maximum AT/T signal together with 
the exponent of a power-law for the relaxation. 



~180 ps. At ripe below 10"'^^ cm~^, the relaxation is still 
power-law (Fig. ^ inset) with an exponent leveling off at 
^-0.5 below 7ipe ~ 1.5 x 10^^ cm~^ that persists out to 
at least ~500 picoseconcP^. 

The plateau region in AT/T is associated with the 
decay of charge into the impurity band into quasi- 
equilbrium. All subsequent relaxation concerns the rear- 
rangement of charge within this band of localized states. 



Relaxation occurs within localized states as the form of 
the THz conductivity does not change after photoexcita- 
tion, as shown by the essentially identical shape of E-field 
transients in Fig. [T^. Therefore tr— >-0 as w— >-0 just as in 
the equilibrium case and as one expects for AC transport 
within localized states. Therefore ct— >0 as w— >0 just as in 
the equilibrium case and as one expects for AC transport 
within localized states. Note also that despite excited 
charge densities of at most 23% of n, the overall scale 
of the transmission change (Fig. [2] inset) never reaches 
more than 0.1% of the unpumped transmission even at 
the highest ripe. This gives additional evidence that the 
decay should be associated with relaxation within a lo- 
calized band of low intrinsic conductivity. 

The strong ripe dependence of the instantaneous de- 
cay rate and power-law is unusual and evidence for 
multi-particle relaxation. Still, even in many systems 
where two-particle recombination kinetics are important, 
the Rothwarf- Taylor "phonon-bottleneck" effect typi- 
cally sets in, which decreases the ripe dependence^. 
In such a case, the phonons created via electron-hole 
pair decay are available to make new electron-hole pairs, 
which limits the net energy flow out of the electronic sys- 
tem and diminishes the role of two-particle reactions as 
a rate limiting step. Effects of this kind apparently are 
not dominant here. The npe dependence is also different 
than relaxation from above gap excitations in other dis- 
ordered semiconductors like microcrystal silicorpSl where 
the decay rate decreases with increasing fipe- There, 
with increasing ripe, the available trapping sites are filled 
up and only slow conduction band - valence band recom- 
bination is available as a decay channel. Obviously, the 
functional dependencies here are very different. 

We propose that the power-law decay and the strong 
tendency towards a smaller exponent at low npe is con- 
sistent with glassy many-body relaxation. After initial 
photoexcitation charges fall back to the impurity band 
into a quasi-equilibrium configuration, which necessar- 
ily occupies higher energy levels where localization is 
weaker. The system is constrained through large poten- 
tial barriers from reaching lower energy configurations 
expeditiously. It is only through rare and intrinsically 
slow many-particle rearrangements that electrons find 
lower energy configurations thereby reducing the con- 
ductivity. These relaxations may occur through Auger- 
like rearrangements^^ of a discrete number of electrons 
or through a coupling to collective electronic modes^®. 

Although we find power-law and not logarithmic decay 
as in the DC experiments on In2 0;^i^, note that power- 
law t"' decay with a <C 1 has a dependence very close to 
logarithmic. Like the logarithmic dependence the power- 
law presumably derives also from an average over a broad 
distribution of relaxation times. Although a power-law 
decay is by definition scale invariant, it is in the sense of 
the small power-law exponent that the relaxation can be 
considered 'slow' and glassy. 

In the present case, despite the decay's scale-free na- 
ture we can compare the rough time scales in our experi- 




FIG. 3: (color) Time-resolved AT/T photoinduced transmis- 
sion dynamics of Si:P sample at different temperatures for 
npe — 1.6 X 10^^ cm~^ and an excitation energy of 0.82 eV. 
The data are normalized to maximum signal. Power-law fits 
are shown as black dashed lines. The inset shows the temper- 
ature dependence of the maximum AT/T signal together with 
the time it takes to decay to 1/e for npe — 1.6 x 10^'' cm~^ 
and npe ^ 6.4 x 10^® cm'^ 



mental range to other natural relaxational scales. At 5 K 
and npe ~ 3.2 x 10^^ cm'^, it takes «200 ps for AT/T 
to decay to 30% of its maximum. This is much longer 
than for instance, the natural scale for phonon assisted 
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relaxation, which should be of order -!-e*^ t •* 
cosecond in this sample at 5 KP. Here cjd is the Debye 

frequency of silicon and Tes is the characteristic Efros- 

1 ^ 
Shklovskii temperature -j^^^- It is also interesting to 

compare to the 'Maxwell time' tm, which is the classical 
expectation for spatial charge relaxatiorP. tm calculated 
from the 5 K DC resistivity^SI p^^ is approximately 0.9 
millisecond {tm = p£)C'/47reieo), which is far in excess 
of any scale we observe. However, it may be possible to 
define an 'AC Maxwell time' from the equilibrium THz 
range resistivity associated with short length scale relax- 
ation. Such a time calculated from older data at our 
central THz frequency (600 GHz)P is 14 ps. The dif- 
ference between these time scales and the experimental 
ones may also be evidence for glassy effects. 

The temperature dependence supports our claim of 
glassy relaxation. Fig. [3] shows the normalized time- 
resolved photoinduced AT/T at various temperatures for 
ripe — 1-6 X 10^^ cm~^. It features a ~1 picosecond rise, 
followed by a nearly constant signal over ~7 picosecond 
after which it displays a temperature dependent power- 
law decay. At the lowest temperature the power-law ex- 
ponent for this ripe is ~-0.6, and becomes about ~-2.5 
above 25 K. The inset shows the dependence of the trans- 
mission maximum for both ripe — 1-6 x 10^^ cm~'^ and 
ripe — 6.4 x 10^® cm^'^. In Fig. Wk we show a compila- 
tion of the temperature dependent fits of the power-law 
decay at different npe. Above a threshold temperature 
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FIG. 4: (color) (a) Exponent of power-law decay vs. tem- 
perature at difTerent photoexcited carrier densities. A low 
temperature saturation region is seen which is temperature 
independent below some npe. (b) Initial decay rate after ini- 
tial excitation vs. npe. The red curve is a power-law fit to 
the intermediate npe region with exponent ~ 1.8. 



of around ^-^25 K and n ~ 1.6 x 10^^ cm~'^, the dynamics 
shows no significant changes with ripe as is expected for 
such a system out of the glassy regime. Below ^25 K, the 
power-law relaxation becomes slower. As noted previ- 
ously the very lowest measured npe exhibits a power-law 
very close to 0.5 at the lowest temperature. Note that all 
measured curves show a low temperature region where 
the exponent saturates. This region shrinks slightly with 
decreasing temperature, but remains finite even in the 
limit of low ripe. This is consistent with quantum tunnel- 
ing as means to relaxation. In contrast, if the relaxation 
was purely due to thermal effects or dependent on the 
energy deposited from the pump pulse, one would expect 
that the temperature threshold for saturation and the 



exponent to go to zero continuously at low T and npe. 

We have shown evidence for a highly temperature and 
npe dependent relaxation consistent with glassy dynam- 
ics. What can we say more precisely about the na- 
ture of the relaxation? In materials where two-particle 
recombination is the principle relaxation channel, one 
expects a — /3npe term in the rate equations and an 
instantaneous decay rate that goes as /3npe. This is 
what has been found in systems like some superconduc- 
tors where electron-hole recombination dominate d^^ l ^^ l 
In glassy systems, simultaneous multi-particle hoppings 
should contribute terms to the rate equations as a high 
order in npe and hence the instantaneous rates would be 
proportional to npe to a power larger than unity. One 
may also expect a super-linear dependence if relaxation 
is driven by excited collective electronic modeS^^l 

In Fig. |4]d we plot the instantaneous decay rate from 
a 10 picosecond time span after the transmission plateau 
vs. npe- Strictly exponential decay would be linear and 
two-particle relaxation quadratic on this plot. The data 
shows a dependence faster than linear and at interme- 
diate npe's goes like the « 1.8 power. This is consis- 
tent with multi-particle relaxation. Due to the previ- 
ously mentioned temperature and npe independence at 
low npe, for a full description it seems essential to in- 
corporate quantum relaxation effects. These are not nor- 
mally considered in relaxation models and theoretical in- 
put is essential for a more rigorous interpretation. 
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